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We study the second quantized -or guiding center- description of the torus Laughlin state. Our 
main focus is the change of the guiding center degrees of freedom with the torus geometry, which 
we show to be generated by a two-body operator. We demonstrate that this operator can be 
used to evolve the full torus Laughlin state at given modular parameter r from its simple (Slater- 
determinant) thin torus limit, thus giving rise to a new presentation of the torus Laughlin state in 
terms of its "root partition" and an exponential of a two-body operator. This operator therefore 
generates in particular the adiabatic evolution between Laughlin states on regular tori and the 
quasi-one-dimensional thin torus limit. We make contact with the recently introduced notion of 
a "Hall viscosity" for fractional quantum Hall states, to which our two-body operator is naturally 
related, and which serves as a demonstration of our method to generate the Laughlin state on the 
torus. 

PACS numbers: PACS 



I. INTRODUCTION 

The discovery of the fractional quantum Hall effeclP 
has led to a series of remarkable theoretical devel- 
opments. Much insight has flowed from principles 
that govern the construction of certain "special wave 
functions"!^!] and their relation with conformal field the- 
ory (CFT)! 3 ^ The credibility of this approach is greatly 
enhanced by the construction of parent Hamiltonians for 
such wave fun ctions, w hich has been possible in many in- 
teresting cases HfflgSSl This very particular class of quasi- 
solvable 11 Hamiltonians consists of Landau level pro- 
jected ultra-local interactions, which enforce the ana- 
lytic properties that uniquely characterize the respec- 
tive ground state. The prime example for such a parent 
Hamiltonian is given by the Vi-pseudopotentialj - which is 
a pairwise (two-particle) projection operator onto states 
of relative angular momentum 1 within the lowest Lan- 
dau level (LLL). Its unique ground state at filling factor 
v = 1/3 is the Laughlin state corresponding to this filling. 

Due to the Landau level projection, the pseudo- 
potential Hamiltonian acts only on the "guiding center" 
degrees of freedom, which exhaust the large degeneracy 
within a given Landau level, and commute with the gen- 
erators of inter-level transitions. (The latter are related 
to the kinetic momenta of the particles, see below. For a 
review of physics in a magnetic field, the reader is referred 
to Ref. [T^l) It is therefore beneficial to make the ac- 
tion of the Hamiltonian on guiding center variables man- 
ifest. This is in particular the case when the Hamiltonian 
is expressed using creation/annihilation operators for a 
set of eigenstates, say, of one of the two non-commuting 
guiding center components, which form a basis for the 
LLL. In numerics, such a second quantized "guiding cen- 
ter" description of the Hamiltonian is essential to make 
use of the reduced Hilbert space dimensionality owing 
to the LLL projection. We illustrate this procedure for 
the cylinder geometry, for reasons that will soon become 



apparent. To this end, we introduce a set of LLL basis 
states as described above, given by: 

4> n (z)=ee-* x2 e-i n2 / r2 (1) 

where £ = e z l r is an analytic function of z = x + iy that 
satisfies periodic boundary conditions in y, appropriate 
for a cylinder of perimeter 2irr (using Landau gauge, A = 
(0, x)). These orbitals are eigenstates of the x-component 
of the guiding center with eigenvalues n/r, where, for the 
time being, we set the magnetic length lg equal to 1. The 
1/3-Laughlin state on the cylinder is then expressed as^ 

iP 1/3 ( Zl ...z N ) = JJ& - 0) 3 xe-^*«I. (2) 
i<j 

With respect to the basis Eq. ([IJ, the V\ pseudo- 
potential takes on the following second quantized form 
(cf, e.g., Ref. HI: 

R (3) 
Qr = ^xcxp(-x 2 /r 2 ) c R - x c R+x 

X 

In the first line, the sum goes over both integer and half- 
odd integer values of R, whereas in the second it goes 
over integer (half-odd integer) if R is integer (half-odd 
integer), such that labels i? ± x are then always integer. 

The one parameter family of models (J3j) share many 
features with one-dimensional (ID) lattice models that 
arise elsewhere in solid state physics, such as transla- 
tional invariance and short ranged (exponentially decay- 
ing) interactions. It is thus not surprising that it has 
recently been proposed to be of use in the absence of 
(proper) Landau level physics, e.g., in flat band solids 
both witfPland withoutP non-zero Chern numbers, and 
in quite general terms in Ref. 1141 

Despite the usefulness of the second quantized descrip- 
tion j3| of the pseudo-potential, it would be very difficult 
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to solve for the zero energy eigenstates of the model in 
this form, or to even know analytically that such zero 
energy eigenstates exist. For this we rely on the original 
first-quantized definition of the pseudo-potential V\ , and 
on the explicitly known analytic form of the Laughlin 
state, Eq. in terms of ordinary position variables. It 
would be highly non-trivial, however, to come up with 
such a first quantized language for the problem Eq. ^ if 
its connection to LLL orbitals were not a priori known. 
This is so because this language becomes available only 
after proper embedding of the degrees of freedom associ- 
ated with the operators c„, 4 in Eq. ^ into a larger 
Hilbert space. In Eq. ([3|, no information is retained 
about the kinetic momenta that determine the structure 
of the Landau levels. Indeed, as Haldane has recently 
shownp^ by making these kinetic momenta subject to a 
different metric from that entering the interactions, one 
obtains a different way to naturally embed the problem 
|3]) into the larger Hilbert space of square integrable func- 
tions. In this setup, Eq. ^ remains unaltered, but the 
resulting wave function loses the analyt ic pro perties of 
Eq. ([2]) that make the problem tractable ! 17 * 18 ! Moreover, 
the solid state applications mentioned initially represent 
yet another way to embed the problem ([3| into a larger 
Hilbert space. 

These considerations show that "interaction only" 
models such as ^ , especially ones that share the "center- 
of-mass conserving" property,^ may enjoy a consider- 
able range applications, but at the same time, may be 
quite hard to solve in general.^ This is chiefly due to 
the fact that the Laughlin state, in its second quan- 
tized/guiding center presentation, is quite a bit more 
complicated than in its analytic first quantized form 
Eq. |2]). While no closed form seems to be known for 
the amplitudes (0\c ni . . . c nN 1^1/3)1 much progress has re- 
cently been made in understanding their structure for the 
cylinder geometry, and for any other geometry in which 
the analytic part of Laughlin's wave function is given 
by a polynomial. Indeed, for Laughlin states and many 
other quantum Hall trial wave functions, these polyno- 
mials have been identi fied a s Jack polynomials, multi- 
plied by Jastrow factors {2212U This allows the amplitudes 
(0|c ni . . . c nN \ipi/3) to be determined recursively. For the 
cylinder Laughlin state, this can be sketched as follows. 
We consider the expansion of Eq. ^ into monomials, 



known recursively, starting from the coefficient of the 



{n k } k 



(4) 



The product in the above equation can be interpreted as a 
state with definite single particle occupation numbers, up 
to a normalization. (The C{„ fc } have the proper (anti- 
symmetry to allow (anti)-symmetrization of the prod- 
uct.) This normalization is readily read from Eq. ([!]). 
We thus hav«Pl 



1/3/r 



{n k } 



c 



{nk} 



• 4jo). (5) 



"root configuration" cj, 



. c t JO) = 11001001 0010. ..) 



through a process known as "inward squeezing" l 2 ° l 21 l 

A remarkable aspect of Eq. ^ is that the dependence 
on geometry, in this case the cylinder radius r, comes 
in only through the trivial normalization factor. This 
is matched by a similarly trivial r-dependence of the in- 
teraction V\. It is quite easy to see that the condition 
that the Hamiltonian Eq. ( 3]) has a zero energy eigenstate 
(which, by positive semi-definiteness, must be a ground 
state), which reduces to Qr^i/^t = Vi?, yields an r- 
independent condition on the coefficients Ci nk \. In this 
way it becomes manifest that regardless of the value or 
r, one is always solving the same problem, which is intu- 
itively clear from the simple analytic form of the Laughlin 
wave function ^ and its trivial r-dependence. It should 
also be emphasized that the the simple r-dependence of 
Eq. ^ is not particular to the Laughlin state. It is a 
direct consequence of the polynomial form of the wave 
function, and carries over without change to any quan- 
tum Hall trial state on the cylinder. 

The situation is rather different for the torus geometry. 
The main purpose of this work will be to get a handle 
on the guiding center presentation of the torus Laughlin 
states. In the remainder of this introduction, we review 
some well known facts that make life more complicated 
on the torus. 

In first quantized language, we pass to the torus by 
introducing periodic boundary conditions in the complex 
plane along two fundamental periods L\ and L2, where 
Li is taken to be real, and ImI/2 > (Fig. [T]). The ge- 
ometry of the torus can be parameterized by r = L 2 /Li, 
the modular parameter. The Laughlin state at general 
filling factor 1 jq then becomes^ 



^[/q{ z l--- z N) =CXp 



i<3 



Li 



(6) 

Here, 8i(z,t) is the odd Jacobi theta-function, and for 
the factor depending on the "center of mass" Z = z\ + 
. . . z jv , which also depends on an additional label I = 
. . . q — 1 corresponding to a choice of basis in the g-fold 
degenerate^ ground state space, we adopt the convention 
of Ref. M 



F e (Z,r) = 



l_ L-q 
2 



(qZ/Li,qr) 



(7) 



The monomial coefficients do not depend on r, and are 



Here, 6*^] (z, r) is the Jacobi theta function of character- 
istics a and b, and L = 2-kL^\yi\Li is the number of flux 
quanta penetrating the surface of the torus. 

Thus, while the Laughlin state is still of the general 
form of a Gaussian factor multiplying an analytic func- 
tion in the complex particle coordinates Zi, the latter is 
not of polynomial form. As a result, to the best of our 
knowledge, there is currently no detailed understanding 
of the structure of the guiding center description of this 
state. By this we mean a general understanding of the 
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FIG. 1. Fundamental domian for torus wave functions, 
coefficients of the analog of Eq. ([5| : 

\rt/ 3 )T = Y, c U}(T)ci N ...ci 1 \o). (8) 

{n k } 

In particular, the r-dependence of the coefficients 
C| nfc j(r) is not of a simple form reminiscent of the r- 
dependence explicit in Eq. Moreover, intuitively, one 
would still expect that these coefficients can be gener- 
ated from the dominance pattern, i.e., 100100100 ... at 
v = 1/3. Indeed, this configuration is still dominant on 
the torus, in the sense that it is the configuration that 
dominates in the thin torus limhP^ (see also Refs. l2l)land 
I26|) . The success of the thin torus approach in determin- 
ing physical properties, such as Abelian and non-Abelian 
statistical^ anc [ presence of gapless excitations^, 
suggests that even on the torus these patterns allow for a 
reconstruction of the full many-body wave function. On 
the other hand, there is no notion of "inward" squeezing 
on the torus, due to periodic boundary conditions. The 
main result of this paper will be the development of a 
machinery for the above mentioned reconstruction of the 
full torus Laughlin state in the guiding center descrip- 
tion, from the thin torus state. Since as an additional 
complication, such machinery can be expected to depend 
non-trivially on r, we first focus our attention on the de- 
pendence of the coefficients C{„ fe }( r ) on the geometric 
parameter. 

As a final remark, we point oufl^that the torus Laugh- 
lin states at v = 1/3 are still the unique ground state 
of the V\ pseudo-potential. Its second quantized form 
agrees with a straightforward periodization of the model 
@, with 



Qh= £ 



E, r x ,2-KlT , > 

;.r • ml.) exp - ' 1 ' '~ 

0<x<L/2 m gZ 



(x+mL) 2 } c R - x c R+x 
(9) 

One sees that for r = z|i2 |/27rr, L2 = iL/r, this reduces 
to the cylinder form Eq. (|3| for L — > 00, and respects 
the periodic boundary condition c n = c n+ L otherwise. 
(Eq. ([9]) is valid for general complex r, though). One 
therefore passes from Eq. ([3| to Eq. ([9| (with imagi- 
nary r) through straightforward introduction of periodic 
boundary conditions (PBCs). Yet the solution of Eq. |9| 
is arguably much less under control. The introduction of 
PBCs is a standard and very useful tool throughout solid 
state physics. We thus expect that a better understand- 
ing of the guiding center description of the torus Laugh- 



lin state will also benefit the solid state application ;] 15 * 16 ! 
mentioned initially. 

The remainder of the paper is organized as follows. In 
section [TT] we construct a two-body operators that gen- 
erated the changes in the guiding center variables of the 
torus Laughlin state with modular parameter r. Sec- 
tions |II A| and |IIB| highlight further formal similarities 
and differences between the cylinder and the torus. Sec. 
II C| presents the heat equation for the r-derivative of 
the analytic Laughlin state. Sec. |IID| introduces a 2D 
to ID mapping, which is our device for embedding low- 
est Landau levels at different modular parameter r into 
the same larger Hilbert space. In Se c. |IIE we derive the 



generator mentioned above. In Sec. |IIF we symmetrize 



this operator and present a byproduct of this study, a 
hitherto unknown class of two-body operators that anni- 
hilate the torus Laughlin state. In Sec. II G we postu- 
late a presentation of the torus Laughlin state in terms 
of its thin torus, or "dominance" pattern, and the class 
of two-body operators generating changes in geometry. 
In Sec. [Hi] we demonstrate the postulate of Sec. II G 
numerically, and work out the relation of our generator 
with the Hall viscosity^!, which we calculate numerically 
as a demonstration of analytical results, comparing the 
resulting data to earlier numerical studies. We discuss 
our results in Sec. |IV[ and conclude in Sec. [V] A small 
Appendix discusses a minor technical detail. 



II. CONSTRUCTION OF THE 2-BODY 
OPERATOR 

A. A final look at the cylinder case 

As motivated above, to establish a machinery that gen- 
erates the full guiding center description of the Laughlin 
state from the root configuration, a natural starting point 
is to get under control how this description changes with 
the geometric parameter t. To this end, we will seek to 
construct an operator that generates changes of the guid- 
ing center degrees of freedom to first order in dr. The 
similar problem for the cylinder, where r is the geomet- 
ric parameter, is comparatively trivial and was already 
addressed in the introduction. For later reference, it is 
instructive to first cast these results in terms of a genera- 
tor of infinitesimal changes in the parameter r~ 2 . Eq. Q 
can be written as 



l^l/3)r' 



= e (r ' ~ r ~ 



)G„-2 



1^1/3)7 



where 



Gr-2 = 2 XX4 C « 



(10) 



(11) 



is the generator of changes in the geometric parameter 
r~ 2 . Note that it is independent of r. We emphasize 



again that ( 10 1, ( 11 ) are very general, and apply to other 
quantum Hall trial states on the cylinder as well. In writ- 



ing (10), we leave it understood that the exponentiated 
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operator generates the change of the guiding center de- 
grees of freedom only; it does not generate in any way the 
change of the LLL orbitals themselves as a function of r, 
Eq. ([I]). We are only concerned here with the change 
in guiding center degrees of freedom, since the object 
of study is the second quantized Hamiltonian Eq. Q, 
in which degrees of freedom associated with kinetic mo- 
menta are not retained. We will thus carefully distin- 
guish from now on between the Laughlin wave function 
tpi/3 = "01/3(^1 ■ ■ .zn]v), which lives in the full Hilbert 
space of square integrable functions over some domain, 
and the ket \ipi/^)r, which lives in an abstract Hilbert 
space denoted C that is isomorphic to the LLL for any 
given value of cylinder radius r. Similar conventions will 
be used below for the torus. In £, therefore, all those or- 
bitals with the same guiding center quantum number n 
become identified, which originally belonged to different 
LLLs corresponding to different values of the parameter 

r [33] 

We note that a similarly universal operator that gen- 
erates changes of the guiding center degrees of freedom 
in respon se to a change in geometry can be obtained on 
the planeP^^ Here, since a geometric deformation by 
means of uniform strain does not affect boundary condi- 
tions, such a deformation is implemented by a change in 
the metric, and unlike in Eq. (10), the operation imple- 



menting this deformation is unitary. 

On the other hand, it is worth pointing out that in 



(10), the lack of unitarity leads to a breakdown of the 
equation in the "thin cylinder" limit r — > 0, in which 
1^1/3) r approaches |100100100 . . . ). The equation re- 
mains valid for arbitrarily small but finite r, where the 
limiting state 1 100100100 ... ) receives arbitrarily small 
corrections, which are, however, important and may not 
be dropped, since they become large under the non- 
unitary evolution facilitated by the exponential opera- 
tor. This is immediately clear from the fact that the 
thin cylinder state is an eigenstate of the one-body oper- 
ator in the exponent. This operator is thus not capable 
of generating the off-diagonal matrix element needed to 
"squeeze" the full many-body wave function out of the 
thin cylinder state, i.e., the root configuration. Eq. ( [To| ) 
is thus not a tool to generate the full cylinder Laughlin 
state out of the root configuration. For the cylinder, how- 
ever, other such too ls are already available, as mentioned 
in the Introduction ! 20 1 21 1 



General considerations for the generator on the 
torus 



We desire to construct an operator analogous to G T -i 
for the torus Laughlin state, which generates changes in 
the guiding center variables of the state in r. This oper- 
ator is thus defined by the following equation: 



Here, G T denotes the operator valued two-component 
object {G Tx ,G Ty ), and V r = {d Tlc ,d Ty ), Note that we 
require that G r is independent of the label I distinguish- 
ing the q degenerate Laughlin states |i/>^(t)), at given 
filling factor l/q and given r. 

To highlight considerable differences with the similar 
problem on the cylinder, we now show that it follows 
easily from these assumptions that, unlike for the cylin- 
der, the components of G r cannot be one-body opera- 
tors. For, if G Tx were one-body operators, we could 
symmetrize each with respect to the magnetic translation 
group. After symmetrization, G Tx would still satisfy 
Eq. fl2| ). This follows from the observation that G T was 
assumed to be independent of £, and that the Laughlin 
states \tp{/ q ( T )) are closed under magnetic translations. 
However, the only one-body operator that is invariant 
under magnetic translations is, up to constants, the par- 
ticle number operator N. Since the \i^x/ q (T)) are eigen- 

states of N, it is clear that no such operators could satisfy 
Eq. @. 

In the following, we will, however, show that G Tx can 
be a two-body operator. 



C. Heat equation for the torus Laughlin state 

We begin by deriving a differential equation for the r 
evolution of the analytic Laughlin wave function Eq. (|6| . 
We have 

d r i>[ /q = e-^« y H(d T F e )f rel +F e d T f rel ) 

where f re i denotes the theta-function Jastrow factor in 
Eq. (|6j) and d T = \(d Tx —id Ty ). The center-of-mass factor 
in the form Eq. Wh is also given by a theta function. 
Independent of £, it satisfies the "heat equation" 

d T F e (Z,r) = J-d 2 z F e (Z,r) = J-d x F e (Z,r) (13) 



4iriq 



Airiq 



with X = Re Z. Since d x leaves the relative part invari- 
ant, the operator (4iriq)~ 1 d x acting on the torus Laugh- 
lin state produces just the first term above in d T ipi/ q . 
The latter can thus be expressed as 



1 



Airiq 



i<3 



d T 6i(zi 



3 1 



i(zi - zj,t) 



V4 



ii 



V T |^ / ,(r)) = G T |^ /fl (r)) 



(12) 



(14) 

It is pleasing that the differential operator on the right 
hand side of the above equation has the form of a two- 
body operator. The are, however, two remaining obsta- 
cles before we can express the change of guiding center 
variables in terms of a two-body operator derived from 
the above equation. First, as defined thus far, the Laugh- 
lin states Eq. ^ for different parameter r do not live in 
the same Hilbert space. In particular, for fixed r the state 
^ is usually viewed as a member of the Hilbert space 
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of square integrable functions over the fundamental do- 
main in Fig. [T] In order to view the differential operator 
in Eq. |i"4| ) as an operator in some Hilbert space, we must 
therefore first embed all Laughlin states for different r, in 
fact all the corresponding lowest Landau levels, into the 
same Hilbert space, since our differential operator can be 
viewed as connecting states with infinitesimally different 
t. The second obstacle is that even with such embed- 
ding, the lowest Landau level will depend on r, i.e., will 
correspond to a different subspace of the larger Hilbert 
space % (to be defined below) for different r. The differ- 
ential operator in Eq. ( 14 ) therefore not only describes 



the change of guiding center degrees of freedom with r, 
it also describes the change of the Landau level itself, 
which we are not interested in. We will therefore find it 



necessary to extract the piece of Eq. ( 14 ) that acts on 
guiding centers only. 



D. Mapping the problem to ID 

We will first address the more technical problem, which 
is the embedding of the torus Landau levels for differ- 
ent t, denoted by C T in the following, into the same 
larger Hilbert space T~L. One natural approach that has 
been emphasized in the recent literature^! is to choose 
an equivalent way to formulate the problem, where the 
fundamental domain remains unchanged and instead the 
metric is deformed. We will return to this point of view 
in Sec. |III[ where we make connection with the Hall vis- 
cosity. 

Here we will choose a different approach, which 
is rooted in the intuition that Landau-level-projected 
physics is effectively one-dimensional. One manifesta- 
tion of this is the form of the "ID lattice" Hamiltonian 
Eq. ([3| that governs the guiding center degrees of free- 
dom. Another is the fact that wave functions in the LLL 
are entirely determined by holomorphic functions satis- 
fying certain boundary conditions. As is well known, 
the values of such functions in the entire complex plane 
are already determined by those on (any interval on) the 
real axis. For this reason we may restrict our study of the 
Laughlin states ^ to the real axis without any loss of in- 
formation. Also, we find it convenient to choose L\ = 1, 
L 2 = r as the fundamental domain for the original two- 
dimensional (2D) wave functions. With this, after re- 
striction to the real axis, all states ^ become elements 
of H = L 2 [0, 1] of square integrable functions within 
the interval [0, 1]. We note that with these conventions, 
the area of the fundamental domain is not preserved as 
we change r. Therefore, we must accommodate for this 
by changing the magnetic length accordingly, such that 
Imr = 2irLl%. This, however, results only in the following 
trivial modification of the wave functions ([6|, 



exp(- 



exp(- 



2 ^ 

k 



following. Clearly, when Eq. ( 14 ) is now restricted to 
y = 0, the operator on the right hand side is a well defined 
differential operator within the Hilbert space % (in the 
usual sense that its domain is dense in %.) 

A preferred basis for the LLL at given r, both within 
the original 2D as well as the ID Hilbert space, is given 
by the following wave functions, 



Xn(z) 



2L\ 



n/L 




(Lz,Lt) 



(15) 



These are eigenstates of the operator exp(^7r a ), where 
TT y is the guiding center y-components. For any r, the 
restriction of these orbitals to the real axis spans a dif- 
ferent subspace C T of the ID Hilbert space H, which is in 
one-to-one correspondence with the lowest Landau level 
at t. 

To see why the orbitals Xn are a natural choice of basis 
in the present context, we observe that the mapping to 
the ID Hilbert space T-L introduces a new scalar product 
between wave functions, defined as usual by integration 
over [0, 1] (instead of integration over the fundamental 
domain in 2D). Eq. (15 1 as written is normalized inde- 



pendent of n with respect to the 2D scalar product, but 
not with respect to the scalar product of fi. However, 
these orbitals are orthogonal in both cases thanks to triv- 
ial considerations of properties under translation in x, 
which are unaffected by the ID mapping. The fact that 
the basis Eq. (151 remains orthogonal, and in particu- 



lar linearly independent, after restriction to the real axis 
makes it manifest that the mapping between the original 
lowest Landau level and its image C T in the ID Hilbert 
space is one-to-one. 

We note that working with y-guiding-center eigen- 
states instead of x (as in our initial discussion for the 
cylinder) leaves the second quantized Hamiltonian invari- 
ant, except for the trivial replacement t — > — 1/t associ- 
ated with the "modular S transformation" . This is due 
to the "S-duality" of the physics on the torus (see, e.g., 
Ref . I55|) . The torus Hamiltonian ([9| was already written 
with reference to the orbitals (15). 



E. 



Definition of a 2-body operator generating the 
deformation of guiding center variables 



We first explain how to relate a result obtained within 
the ID framework introduced above to the desired one, 
which uses ordinary conventions based on a Hilbert space 
equipped with the standard 2D scalar product. 

Suppose we have an operator G Tx {G Ty ) that generates 
the change with t x (r y ) in the coefficients C{ nfc }(r) in the 
expansion of the Laughlin state, 



which is inconsequential since we work at y = in the 



• Xn N (zN,r), 
(16) 
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where X„(z,r) = A/" n (r)x n (^,r), 7V„(r) being the fac- 
tor that normalizes the state \ n with respect to the 
ID scalar product, i(</>|V')i — Jq dxcj)*(x)4>(x), i.e., 

1/2 

■Mi = i{Xn\Xn)i , and we will often leave the r- 
dependence understood. Likewise, we have dropped the 
label I for now, which is just a spectator in the "heat 
equation" (14). A denotes anti-symmetrization in the 
indices rik- Eq. (16) is a member of the subspace C T of 
% as defined in the preceding section. We may now map 
the state ( 16 ) to the abstract Landau level Hilbcrt space 
£ as discussed in Sec. |II A| by applying a projector which 
"forgets" the degrees of freedom associated with kinetic 
momenta. This situation is represented by the diagram 
in Fig. [2] If we perform this projection orthogonally with 



t 



J flow of Eq. (14) 
L-r > Cn 



It 



h' 

4->jv„(/)4 



p e s: g t 



G T £ 



FIG. 2. Commuting diagram displaying the various Hilbert 
spaces and sub-spaces denned in the main text, and operators 
acting between them. The top segment shows lowest Landau 
levels £ T , C T ' at different modular parameter that, using the 
2D to ID mapping defined in the text, have been embedded 
into the same larger ID Hilbert space H. At the same time, 
each Landau level is isomorphic (through embeddings I T , I T i) 
to the same finite dimensional "abstract" Landau level space 
£, in which only the guiding center degrees of freedom are rep- 
resented. The generator G T of changes in the guiding center 
degrees of freedom with r is first constructed using the nor- 
malization conventions of the ID Hilbert space. It is related 
by a similarity transformation to the operator G T , which gen- 
erates the analogous changes for the normalization convention 
of the usual 2D Hilbert space. In the horizontal direction, we 
have mappings between states defined for values of the mod- 
ular parameter. The upper line is defined through the flow of 
Eq. (14), which describes precisely the change of the Laugh- 



lin state, restricted to the real axis. The lower line represents 
the corresponding change in guiding center degrees of free- 
dom, given by Eq. (41 1. The operator G T is constructed such 



that the diagram commutes. 

respect to the ID scalar product, we obtain a ket 

l&/«to> = E C {nk} (r)cl...cl N \0). (17) 

{»*} 

By definition, we then have 

V T $i/,(r)> = G T |&/,(r)) , (18) 



where we assume G T = (G Tx , G Ty ) to be of the form 



mm Tin 



In the end, one wants to do the projection of Eq. ( 16 1 or- 



thogonally with respect to the original 2D scalar product. 
This gives 



hM 1 ")) = E 
{«)=} 



•4J0) (20) 



where C{, lk } = J\f ni . . ■ N'n N Gs nk \ from the change of nor- 
malization, cjj — > Af n c „- This implies the relation 

llMT)) = e ^oln[^(r)]ct Cn| ^ /g(r)) (21) 

From this last line, we obtain that the desired operator 



G r defined by Eq. ( 12 ) is related to Eq. ( 19 ) vi 



,36 



_ V VtM,(t) t 

A/" m (r)Mn' (t , r , 

+ 2^ .r,s.r ^ 



mm'nn' 



mm'nn' clnC nl /C n C n ' . (22) 



With this we have completely relegated the solution of 
the problem to the ID Hilbert space. We point out that 
the ID mapping described above may generally provide 
an efficient way to calculate the matrix elements of opera- 
tors acting within the lowest Landau level on the torus.^ 
In this case, Eq. (22) will apply without the r-derivative 



part. The explicit form of N n (r) will be given below. 

We now define the operator I T which injects the ket 
IV'i/gW) m t° A- £ by sending aJjO) to Xn( T )- Thus 



1pl/g = Irl^l/q) 



(23) 



For the time being, we work at fixed t x . Using the 
heat equation (14) with d T = —id Ty and differentiating 



Eq. (23), one obtains 



dr y lpl/q = iMl/q = {drJrMl/q) + It G Ty \lJj 1/q ) , (24) 

where A denotes the differential operator on the right 
hand side of Eq. (14), and we also used Eq. (18). 

For Re r = Re t\ it is easy to see the P T d T I T = 0, 
where P T is the orthogonal projection operator onto C T 
(we work in the ID Hilbert space now, and will al- 
ways refer to its scalar product when not stated oth- 
erwise). To see this, it is sufficient to observe that 
(Xm\d Ty Xn) — for all m, n. This follows from the fact 
that (Xm{T)\Xn{T')) = <5 m ,„ (XmW \Xm (t')) is always real 
for Rer = Rer'. Thus, acting on the last equation with 
P T , we get 

P T AP T ip 1/q = P T AP T I T \ib 1/q ) = I T G Ty \i> 1/q ) . 

where we have also inserted P T before ipi/ q € C T . Since 
we only care about how the operator G T acts on these q 
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states, for which we have the last equation, we may thus 
define this operator though the identity 



G T =I- l P T &P T I T . 



(25) 



Next, G\ is the con trib ution coming from the differential 
operator d\ in Eq. ( 14 ) . Note that after normal ordering, 
the square of a singlebody operator still contains a single 
body operator. We thus get the following result: 



The last equation expresses that the matrix elements of 
G Ty are just those of the differential operator A restricted 
to the LLL subspace C T . These can thus be calculated 
straightforwardly by evaluating the standard expression 
for two-body operators: 



G 1 = (f) 2 [El^«+ E |^' ! 



CO CO C ni C ™i C n2 C "2-l 



G 



1 ,1 

dx I dx' X* m (x)X* m '( X ') A X*n'(x')X* n (x) ■ 

o Jo 

As a last step, we calculate G Ty by fixing the normal- 
ization convention for single particle orbitals in accor- 
dance with the usual 2D scalar product, as displayed in 
Eq. (22). We may then obtain the generator for changes 



(34) 

Note that S 2 ^ (S^) 2 , owing to the fact that P T d\P T ^ 
{PrdxPr) 2 - Finally, G2 relates to the ^-function part of 
Eq. ( 14 1. Eq. ( 26 ) can be evaluated by expanding the fac- 



tors in the integrand, which are all periodic in x, x' , into 
Fourier series. For the 90/0-terms, this can be done by 
contour integration and using known properties of 9 func- 
tions. Straightforward but tedious calculation allows one 
to express G2 through rapidly converging, albeit multiple 
sums, 



in t x simply by studying the analytic properties of the 
coefficients C{„ fc }(r) in Eq. (20). As shown in Appendix 
|A"| one has 



G 2 



%^F £ 4 i ct 2 c„ 4C „3 + ^iV(iV-l), 



d r x C {n k } = - id r y C {nk }+i 



N 
4% 



(27) 



and we have defined the function 

2tt 



(35) 



We can thus let G T 



-iG Ty +i 



Moreover, Eq. (|27|) 

N/4 



follows from the fact that Gi ni \/T v is holomorphic in 
r. We may use this insight to conveniently redefine the 
normalization of the Laughlin states via 



E ( 

n— 713 — m mod L 



1 - e 



2i7rm ' 



^TLKni+nJ/L+h] 2 ^ e iTvrL{7i 1 /L+l 1 ) 2 y 

h 



1p'l/ q (zi,---, z N,T) 



■,ZN,T) 



(28) 



-i7rrL[tn4+n)/L+^4] 2 ^* g27r rL (77.4 /L+Z4) 



The corresponding generator for changes in r y is then 

given by G' Ty = G Ty — -^r ■ In the following, we will always and the (T _ dependent) constan t 



refer to the normalization convention (28). Dropping all 
primes, we then have 



C 



~iG T = G T . 



(29) 



1/2 A^, 

-1/2 ft 



d'Z — 7T 2 ] 



(36) 



(37) 



With the ket \ipi/ q (r)) now referring to Eq. (28), 
\ 1 Pi/q( T )) i s then holomorphic in r, and we have 



dr\A/ q (r))=G T \^ 1/q (r)). 
We present our final result as 

1 



G T — Gg 



Gi + qG 



2 • 



(30) 



(31) 



Here, the first term corresponds to —i times the one-body 
operator in the y-component of Eq. (22), plus the shift 



of iN / 4-Tj, shown in Eq. ( 27 ) . Defining the functions 



S° = V(27ri[ZL + n])V 



-2irLT y (l+n/L) 2 



(32) 



In the above, the Kroneckcr i5 is understood to be peri- 
odic, enforcing identity n\ + n-i — ns + mod L. 



F. Symmetry considerations 

The operator G T defined in the preceding section is 
manifestly invariant under magnetic translations in the 
x-direction. In the basis we chose here, this is tanta- 
mount to the conservation, modulo L, of the "center- 
of-mass" operator X)n ric n c ™- On ^ ne °th er hand, the 
operator is not invariant under magnetic translations in 
the r-direction, which amounts to an ordinary shift of 
the orbital indices. As already pointed out in Sec. |II B[ 
the symmetrized operator 



the normalization factors defined above correspond to 
(Af n )- 2 = y/2L/r y S°. We thus get 



G T 



1 L ~ l 

71 = 



T T l G T (Tl) n 



(38) 



9rX 



2^ S° 



1 S 2 

" Cn = ~ AiriL ^ C " C ™ ' ^ 



2 
_ 



also satisfies Eq. ( 12 ), where T T generates magnetic trans- 



lations in r. This is a trivial consequence of the fact that 
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the |'!/ , i/ q ( T )) transform among themselves under T T , and 
all satisfy Eq. ( 12 1. Likewise, each term on the right hand 
side of Eq. (38) satisfies Eq. (12). We may thus define 



the L — 1 linearly independent 2-body operators 



D r . 



G T , Bym -T?G T (T}) n , 



0...L-2, 



(39) 



that all annihilate each of the q-fold degenerate Laughlin 
states, 



D n \ip{/ q )=0 



(40) 



We note that the D n are not in any obvious way related 
to the operators Q* r Qr of the pseudo-potential Hamil- 
tonian, with Q_r given by Q. Indeed, the D n have a 
non- vanishing single-body term, whereas the Q r Qr. do 
not. The D n thus represent a new class of two-body op- 
erators that annihilate the torus Laughlin states (in the 
absence of quasi-holes) . For q = 3 and various values of 
particle number N, we have verified that the property 
( 39 1 characterizes the q = 3 Laughlin states uniquely. 

Note that the single-body contribution to G Tj3ym is pro- 
portional to the particle number, as explained in Sec. 
|II B| This term can thus be replaced by a constant when 
acting on the Laughlin state, and hence can be ignored 
altogether in practical calculations, where the real part 
of this constant is usually adjusted to fix the normaliza- 
tion of the state (see below), and the imaginary part only 
affects the phase convention. For the same reason, we do 
not need the value of the r-dependent constant C defined 
in Eq. (37) for the purpose of practical calculations. 



G. 



Presentation of the Laughlin state through its 
thin torus limit 



In the following, we will generally identify G T with the 
symmetrized operator G T . Bym discussed in the preceding 
section, without carrying along the " sym" label. Putting 
the results of Sec. |II E| in integral form, we have, via 
Eq. p|, 



|^ /g (r'))=Pe^' G ^|^/ 9 W), (41) 



where P means path ordering. The integral in Eq. (41) 
should be interpreted as a complex contour integral, 
where the result is independent of the path connecting 
t and t'. This is so since by construction, G T generates 
exactly the change with r of the guiding center coordi- 
nates of the states in Eq. (28), which are single valued 



functions of r. (This requires that we carry along all the 
r-dependent c-number terms mentioned in the preceding 
section.) 

We may also want to add, possibly different, real con- 
stants to G Tx and G T , such that the normalization of 
the Laughlin state is preserved under the evolution with 
these operators. When evaluating Eq. (41) iteratively, 



this simply corresponds to normalizing the state at each 
step. We denote the accordingly modified operators by 



G^: and ,, and introduce the operator valued 1-form 
dG 



N 



G^ dr x + G' v dr v . We may then write 



~<N 



Pe 



\ti /q (r)) N , (42) 



where the subscript N denotes normalized Laughlin 
states. We are now interested in the thin torus limit 
t — > zoo, in which \iP{/ 3 (t))n approaches the ket 
1 100100100...), 14 or one related to the latter through 
repeated action of T T . Here, the labels 100100100. . . are 
occupation numbers in the basis ( |15[ ). Given our earlier 
discussion for the cylinder, it cannot be taken for granted 
that Eq. (42) remains well-defined in this limit. On the 



other hand, it may seem plausible that this is the case, 
since the operators G^ , G^ do generate off-diagonal ma- 
trix elements when acting on the thin torus state, unlike 
the case of the cylinder. It thus seems feasible that the 
full Laughlin states at arbitrary r admit the following 
presentation in terms of their respective thin torus limit, 

\^x/ q {r))N = Pe^ dG ™|...100...100...100 ) , (43) 

where the pattern on the right hand side denotes one of 
the q thin torus patterns at filling factor l/q. The cor- 
rectness of the above assertion remains non-trivial, how- 
ever, as the t' — > oo limit in Eq. ( 42 ) must be taken with 



care. In the next Section we provide numerical evidence 
for q — 3, demonstrating the above relation for various 
particle numbers N. We thus find that the full torus 
Laughlin state may be generated from its given thin torus 
limit via application of the above path-ordered exponen- 
tial involving the two-body operator constructed here. 
We conjecture that this is true for general q. An appli- 
cation demonstrating this technique will be discussed in 
the following. 



III. APPLICATION: HALL VISCOSITY 

As an application of our findings in Sec. [TTJ we use 
Eq. (43) (or the differential form Eq. (30)) to calculate 



the v = 1/3 torus Laughlin state along a contour in 
the complex r-plane, starting from the thin torus limit 
at r = ioo. As a physical motivation for calculating 
the Laughlin state along such contours, we will be ask- 
ing how the Hall viscosity^] evolves along such contours. 
This quantity is naturally related to the main theme of 
of our paper, i.e., changes of the Laughlin state with 
changes in geometry. The notion of a Hall viscosity of 
fractional qua ntum H all liquids has generated much in- 
terest recently! 32 ! 34 * 38 ! expanding earlier worlP^ n inte- 
ger quantum Hall states. In particular, in an insight- 
ful paperpS Read has derived a general relation between 
the viscosity of a quantum Hall fluid and a characteristic 
quantum number s, which can be interpreted as "orbital 
spin per particle" and is related to the conformal field 
theory description of the state in question. Here we only 
give a brief summary of the relevant definitions, following 
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closely Ref. [331 to which we refer the interested reader 
for details. 

We denote the fourth-rank viscosity tensor of the fluid 
by Tjajjcd, where we are interested in the case of two 
spatial dimensions. In a situation with no dissipation, 
only its anti-symmetric or "Hall viscosity" component 

^abld = ~ r lcdib mav be non-zero, and this is possible only 
when time reversal symmetry and the symmetry under 
reflection of space are both broken. This is the situation 
in a magnetic field (where in a constant field, only the 
product the these two symmetries is unbroken). 

We now consider a system with periodic boundary 
conditions defined by two periods Li = (£i,0) and 
L 2 = (LiT x , LiTy), and Hamiltonian 



1 N 



Plll ^2^2V(\\x ij +mL 1 + nL2\\g)P 1 . 



(44) 



Here, ir a is a component of the kinetic momentum, P LLL 
denotes LLL-projection, and we have introduced a met- 
ric g a b- We have also introduced the "periodized" ver- 



sion of a potential V that depends on Xi 



only via 



llajjjllg = g a b%ij%ij with Xij — Xi — Xj. We follow Ref. IM1 
and parametrize the metric via g(X) = A T A, A = exp(A), 
where A can be viewed as a coordinate transformation 
that transforms the identity metric into the metric g. 
Clearly, g is invariant under A — > RA, where R is a rota- 
tion matrix. Since A can be interpreted as being propor- 
tional to an "infinitesimal version" of A, whose rotational 
component is just its anti-symmetric part, we may fix this 
rotational degree of freedom by requiring A to be sym- 
metric. Then, the Ha ll viscos ity of the ground state of 
Eq. ( 44 ) can be related 3 - 4 * 39 * 40 ! to the adiabatic curvature 



on the space of background metrics, here parameterized 
by the symmetric matrix A. Specializing to g = id, we 
have: 



J A] - --TP 
I abed ~ y r ab;cd 

where F ao;c d is the Berry curvature 

Fab-cd = -2lxa{d Xab i)\dx ci i))\ id , 



(45) 



(46) 



and tp denotes the ground state of Eq. (44 1. F a b- C d clearly 



has the anti-symmetry of Vatldi an( ^ ^ ^ s a ^ so symmetric 
in the index pairs ab and cd. Furthermore, at least in 
the thermodynamic limit of large Li, one would expect 
Vabcd t° acquire full rotational symmetry. In two dimen- 
sions, this requires the trace r)^ c to vanish, where we 
use the sum convention, and similarly for the first in- 
dex pair. (In higher dimensions, rotational symmetry 
requires r/" 4 ' to vanish identically). Moreover, in an in- 
compressible fluid, the strain tensor u a (, must be trace- 
less. Therefore, since the viscosity couples to the rate 



of strain u a b via iJabld^cd to give a viscous contribution 

(A) 

to the stress tensor, only the traceless part of ry ab L is of 
interest. It therefore makes sense to restrict our atten- 
tion to traceless A a b, corresponding to volume preserving 
coordinate transformations. Requiring F a i,-cd thus to be 
anti-symmetric, symmetric in the first and second pair, 
as well as traceless, in D — 2 the associated curvature 2- 
form F = \F a i,. c( i dX a b A dX c d can only depend on the fol- 
lowing two independent linear combinations of 1-forms, 
dAn — d\22 and <iAi2 + d\2\- Hence it must be propor- 
tional to their product!^ 



F = 



1 



s (d\u - dA 22 ) A (d\i2 + d\ 2 i) , 



(47) 



and we introduced a proportionality factor — s/2 whose 
physical meaning will be given below. The above expres- 
sion in Eq. ( 45 ) gives 



Vabld = V (A) {Sadtbc + S bc e ad ) 



with 



V 



(A) 



-snh 



(48) 



(49) 



where n — N/V is the particle density, s — s/N, and 
we have restored a factor of h. As shown in Ref. [35J in 
the thermodynamic limit the parameter s is quantized 
and can be identified with the average orbital spin per 
particle, which is related to the conformal dimension of 
the field describing particles in the conformal field the- 
ory description of the state. It is further related to the 
topological shift on the sphere, S, of the underlying state 
via s = S/2. For the Laughlin 1/3 state, s = 3/2. 

We now consider fixed boundary conditions described 
by r, and introduce a metric that corresponds to the 
infinitesimal transformation 



dX 



1 



-dr y dr? 



2t v \ dr x dr y 



(50) 



It is not difficult to see that the corresponding metric 
change is equivalent to changing the modular parameter 

T tO t' = T 

as 

Ns 



dr x + idr y . We may thus rewrite Eq. ( 47 ) 



- ^ dr x A dr y . 



(51) 



To each A can be associated a r', where A = exp(A) is the 
coordinate transformation that changes the r-boundary 
condition into a r'-boundary condition, where 



ALj ■ AL 2 + iKLx x AL 2 



(52) 



For fixed r, we now parameterize A, and thus the metric, 
by t' . (Note that the right hand side of Eq. (52 1 can be 
viewed as a function of A and r.) Eq. (511 then implies 
that 



4r 2 

s = -^-lm(d T ^\d Tly 4)) . 



(53) 
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We emphasize that in the above, j/j = i/j(t, <7t(t')) always 
satisfies the same boundary condition defined by r, and 
depends on t' only through the metric. At the same time, 
-0(t, g T {T r )) is related to iP(t' , id) by the unitary transfor- 
mation Xn(T,g T {r')) -> Xn(T',id), with Xn(T,3( T ')) the 
deformed version of the state ( p~5| ) in the presence of the 
metric g(T r ). However, for fixed r, the ip(r, g T (T r )) live 
in the same Hilbert spaced independent of t'. The ad- 
vantage of introducing both r and r', where the former 
describes boundary conditions, and the latter describes 
the "true geometry" of the system, is that we may restrict 
ourselves to metrics g T (r') m the vicinity of the identity 
(corresponding to t' close to r), such that Eq. (451 is 
directly applicable. 

We now consider -0 = ^>f/ 3) the normalized Laughlin 
1/3 state (where we suppress labels r, r', and I). Wc 
have the expansion 

</3=E c, {na( r ')IK}) 9 , (54) 

{n k } 

and \{nk})g) is short for the Slater determinant 
•A. \Xm(zi,T,g(T')) ■ ... ■ Xn N (zN,T,g(r'))]. We write 
— 21m{d T 'ip\d T 'ip) = V T ' x A where 



1.5r 



A = 



{«fc} 



C? nk} \ 2 9 ({n k }\V T ,\{n k }) g 



(55) 

is the Berry connection. It turns out that in the first 
term, which describes the change of the LLL basis with 
t', 9 ({nfc}| V T / is independent of {n k }, and con- 



tributes a constant 1/2 to Eq. (53). 3D The second term 
depends on the changes of the Cs nk y with r', which we de- 
scribed in the preceding section. We first assume the gen- 
eral situation where this change is described by Eq. ( 12 ) 



with two generators G Tx and G Ty that are not necessarily 
related and that do not necessarily preserve the normal- 
ization of the state. It is straightforward to show that 
the contribution from the second term then leads to the 
following connected expectation value, 



N 



= i [{G\G Tv - Gl y G T J - (Gl)(G Ty ) + (Gl y )(G T J 

(56) 

where expectation values on the right hand side are taken 
in the state Eq. (54 1. The last two terms take care of 



the normalization, and will cancel if both operators are 
anti-Hermitian (describing unitary evolution), in which 
case the expression reduces to the expectation value of a 
commutator. Note also that the expression is invariant 
under constant shifts of any of the two operators. We now 
specialize to the case where these operators are related 
by Eq. ([29]). Plugging Eq. ([55]), Eq. ([56] into Eq. ([53]), 
this gives 



4r 2 

— |AG T 
N 1 



(57) 



Legend Error 

5p 2.5e-5 

7p 3.0e-4 
— 8p 7.3e-4 
9p 1.5e-4 




FIG. 3. Average "orbital spin per particle" s as calculated 
from Eq. (53 1, for the v — 1/3 torus Laughlin state at r 



generated via Eq. (581, using the procedure described in the 
main text and below. We start with the thin torus stat e at 
relatively large but finite r', and iteratively solve Eq. (301 



using the 4th order Runge-Kutta method. Final state errors 
compared with exact diagonalization at r = i are shown for 
5,7, and 8 particles for r' = 30i and step size dr = 0.05i, and 9 
particles for r' = 40i, dr = 0.025i. Errors for different dr and 
t' are quoted in the main text, the plots are indistinguishable. 
Crosses denote the value of s — 1/2 obtained from the exactly 
diagonalized Laughlin state in Eq. (531, for comparison. The 
errors of s are at or smaller than 10 _D in these cases. Errors 
decrease further with larger initial r' and smaller step size. 



where \AG T \ = ((GtG r )-(G]:)(G r )) 1 / 2 is the variance of 
the operator G r in the state ^1/31 an d is manifestly posi- 
tive (the Laughlin state at r certainly being no eigenstate 
of G T for any r). As stated above, for the Laughlin state s 
is expected to approach 3/2 in the thermodynamic limit. 
This has been checked in Ref. |MJ by calculating torus 
Laughlin (and other) states by exact diagonalization of 
parent Hamiltonians, and computing the Berry curva- 
ture by taking overlaps between such states for different 
r (or A). Here we will consider the same problem both 
as a demonstration and a consistency check of the results 
presented in the preceding section. To this end, we cal- 



culate the Laughlin state from the presentation (43), or 



by numerically integrating the differential equation ( 30 ) 



with thin torus initial conditions, and then computing s 
from Eq. (57). Note that both steps of the calculation 
make use of the two-body operator G T . In particular, 
our results will confirm the accuracy of Eq. ( 43 1 , which 
may be written more carefully as 



\iP[,Jt)) n = lim Pe^' dG - |... 100.. .100.. .100 ). 

(58) 

Evaluating the expression on the right for some large but 
finite t' is equivalent to integrating Eq. ( 30 ) (and nor- 



malizing the result), where the thin torus limiting state 
defines the initial condition at t'. This obviously intro- 
duces some error compared to the full Laughlin state at 
the initial value r', hence also at the final value r. Since 
it is not clear a priori how this error behaves in the limit 
of large imaginary r', possible pitfalls are that the limit 
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in Eq. (58) is ill-defined, or that it is well-defined but 
does not agree with the Laughlin stated Our results, 
however, give strong support of Eq. (58). 

Fig. [3] shows the results for the value of s — 1/2 from 
this method for q — 3. Beginning with the thin torus 
state 1 100100100 ... ) at large imaginary t', we evolve 
the state down to r = i, i.e., a torus of aspect ratio 1, 



1.5 



integrating Eq. ( 30 ) using the classical 4th order Runge- 
Kutta method. We normalize the state at each step. For 
particle numbers N = 5 to N = 9, we have observed 
that the error of the state obtained at r = i, compared 
to Laughlin state generated from exact diagonalization of 
the Vi Haldane pseudopotential, becomes systematically 
smaller with increasing initial t' and decreasing step size 
dr. The observed state error at r = i has been on the 
order of 10~ 6 for N = 5, dr = 0.01, and r' = 28, and 
on the order of 10~ 4 for N = 9, dr = 0.025, and r' = 
40. Generally, larger ./V requires larger r' for the same 
accuracy. Fig. [3] shows s - 1/2 for various N using our 
method, whereas crosses denote isolated points for which 
values have been obtained from exact diagonalization for 
comparison. One sees that the expected value of s— 1/2 = 
1 is always approached rather closely for r = i, though 
it deviates from this value for |r| noticeably larger than 
1. The crossover where notable deviations from 1 set in 
is pushed to larger |t| with increased particle number, as 
expected. However, the value of s is found to be much 
more constant, and close to its expected thermodynamic 
limit, when instead of varying the modulus of r we vary 
its phase at |t| = 1, even for five particles, as shown in 
Fig. HI Data were obtained by continued integration of 
Eq. ( 30 ) away from r = i along a contour where r = 
exp(w)7 s — 1/2 remains close to 1 except for angles 
approaching tt. These observations are consistent with 
the exact diagonalization data published in Ref. for 
larger particle number. 



IV. DISCUSSION 

In the preceding sections, we considered the change 
in the guiding center variables, with modular parameter 
r for the torus Laughlin states. Within a given Lan- 
dau level, the guiding center coordinates fully specify 
the state. We have shown that this change is gener- 
ated by a two-body operator G T , which we have explicitly 
constructed. We have demonstrated numerically that by 
means of this two-body operator, the Laughlin state for 
any modular parameter t can be generated from its sim- 
ple thin torus (r = ioo) limit. The ability to generate 
the full torus Laughlin state in this way may be com- 
pared to squeezing rules that follow from the Jack 



nomial structure of this state in other geometries! 20 * 21 ! 
From a practical point of view, however, our method still 
requires integration of a first order differential equation. 
While this requires some compromise between accuracy 
and computational effort, the added benefit is that in the 
process of the calculation, the Laughlin state is generated 
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FIG. 4. Average orbital spin per particle s, calculated from 



Eq. (53 1, for the v = 1/3 torus Laughlin state. The state has 
been evolved out of the thin torus limit first down to r = i 
as described in the caption of Fig. [3] and then to r = e %e 
using the same method. The step sizes used are d8 — O.Olrad 
for 5 particles and d9 — O.OOlrad for 8 particles. The state 
difference with the exactly diagonalized Laughlin state at the 
last step is listed in the figure. Data stay close to the value 
s — 1/2 = 1 expected in the thermodynamic limfp^for a wide 
range of angles 8, as first observed in Ref. |3U from exact 
diagonalization for larger particle numbers. 



along an entire contour in the complex r plane, rather 
than just for a single value of r. It is thus likely that 
whenever a moderate error can be tolerated, but many r 
values are of interest, our method may become competi- 
tive compared to numerical diagonalization. As a demon- 
stration of these features, we have produced results relat- 
ing to the Hall viscosity that are similar to those of Ref. 
|M](and are expected to be identical within numerical ac- 
curacy for identical particle number, which we have not 
yet studied). The Hall viscosity is itself deeply related 
to our main them e of stu dy, i.e., geometric changes in 
the Laughlin stateJ^MMl anc | we have discussed its pre- 
cise relation to the generator constructed here (Eq. (57)), 
following Ref. M 

We note that one key ingredient of our procedure is to 
embed different torus Laughlin states, which are related 
to one another by the application of strain, into the same 
Hilbert space. For this we make use of a dimensional re- 
duction that is made possible by the analytic properties 
of lowest Landau level wave functions on the torus. We 
argued that this mapping may be useful in other con- 
texts. However, recent work on Hall viscositjEH achieves 
the same embedding by a different method, which is to 
introduce a metric describing the effect of strain, rather 
than a change in boundary conditions. We conjecture 
that if we had used this method in Sec. |Hj we would 
have directly obtained the symmetrized version of our 
operator G T . In this way, however, we would not have 
obtained the family of two-body operators given in Sec. 



II F whose members annihilate the torus Laughlin states. 

While primarily, we have been working in a finite di- 
mensional Hilbert space that represents guiding center- 
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coordinates only, the operator defined in Sec. [TT] also 
naturally acts within the full Hilbert space, which can be 
viewed as the tensor product of the degrees of freedom 
for the guiding centers and the dynamical momenta, re- 
spectively. Within this larger, physical, Hilbert space, 
the operator G T generates the change in guiding center 
degrees of freedom associated to a change in the torus 
geometry, but not the corresponding change of the Lan- 
dau level. As pointed out recently by Haldanep^ the 
Laughlin state may be generalized by the introduction of 
a geometric parameter that describes the deformation of 
guiding center variables in response to a change in the 
"interaction metric" . The so deformed Laughlin state is 
still the exact ground state of an appropriately deformed 
Hamiltonian. The operator that we have constructed can 
thus also be viewed as generating the change of the torus 
Laughlin state in response to a change of the interac- 
tion metric, i.e., the change in ground state for the cor- 
responding family of deformed pseudo-potential Hamil- 
tonians. For the disc geometry, this pro blem has been 
addressed from different angles previously ! 18 * 34 ^ 

We conjecture that the observations made here are 
not limited to Laughlin states, but can be generalized 
to other quantum Hall states as well. Indeed, a great 
wealth of model wave functions is obtained from con- 
formal blocks in rational conformal field theoriesP For 
conformal blocks on the torus, the dependence on the 
modular parameter r can be described by Knizhnik- 
Zamolodchikov-Bernard (KZB) type equations.^ We ex- 
pect therefore that our approach can be generalized to 
other trial wave functions related to conformal field theo- 
ries. The details of such generalizations are left for future 
work. 



V. CONCLUSION 

In this work, we have shown that geometric changes 
in the guiding center coordinates of the torus Laughlin 
state are generated by a two-body operator. We have 
demonstrated that the equation that governs the evolu- 
tion of the torus Laughlin state as a function of the mod- 
ular parameter r can be continued into the thin torus 
limit. This gives rise to a new presentation of the torus 
Laughlin state in its second quantized, or guiding center, 
form. This presentation allows one to calculate the torus 
Laughlin states in terms of a simple thin torus or "domi- 
nance" pattern by means of integration of the flow gener- 



ated by the two-body operator defined in this work. This 
operator hence realizes the adiabatic evolution of the sim- 
ple thin torus product state into the full Laughlin state 
on regular tori. To demonstrate this, we have numeri- 
cally compared both the Laughlin state generated from 
this method, as well as the Hall viscosity derived from it, 
to exact diagonalization results. While the demonstra- 
tion of our new presentation of the torus Laughlin state 
rests in part on numerics, we defer more detailed analytic 
studies to future investigation. 
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Appendix A: Analytic properties of coefficients 

For definiteness, we will refer to the Laughlin state 
^ij q using the normalization conventions Q. The 
coefficients Cs nh \(r) defined in Eq. (20) then imply the 
following expansion of the analytic Laughlin state, 



{n k } 



■ Xn N (z N ,T) . 

(Al) 

Here, as before, the symbol A denotes anti- 
symmetrization, and single particle orbitals \ n 
are defined in Eq. (15). We define new orbitals 



Xn( T ) — T y Xn( T ) that are holomorphic in r, as is 
the Laughlin state i>\[ q (j). Hence, by acting with 

df = \{d Tjl + idr y ) on Eq. (Al), we obtain 



°= E [^(cw(t)a 



N/4 



"4x^(zi,t>. . ■■x'„ n (zn,t) . 

(A2) 

The linear independence of the orbitals x' n ( T ) an d of the 
associated many-particle Slater determinants then im- 
plies 



d,(c {nk} (r)/r y N ^ 



0. 



(A3) 



i.e., the quantities C{„ fc } { T )l T y^ are holomorphic in r. 
Eq. (27) follows immediately from Eq. (A3). 
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